Effects of decoherence and imperfections for quantum algorithms 
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Abstract. We study effects of static inter-qubit interactions and random errors in quantum gates on the 
stability of various quantum algorithms including the Grover quantum search algorithm and the quantum 
chaos maps. For the Grover algorithm our numerical and analytical results show existence of regular and 
chaotic phases depending on the static imperfection strength e. The critical border e c between two phases 

drops polynomially with the number of qubits n q as e c ~ n q 3 ^ 2 . In the regular phase (e < e c ) the algorithm 
remains robust against imperfections showing the efficiency gain e c /e for e > 2~ 71q / 2 . In the chaotic phase 
(e > e c ) the algorithm is completely destroyed. The results for the Grover algorithm are compared with 
the imperfection effects for quantum algorithms of quantum chaos maps where the universal law for the 
fidelity decay is given by the Random Matrix Theory (RMT). We also discuss a new gyroscopic quantum 
error correction method which allows to reduce the effect of static imperfections. In spite of this decay 
GYQEC allows to obtain a significant gain in the accuracy of quantum computations. 
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In realistic quantum computations pQ the elementary 
gates are never perfect and therefore it is very impor- 

, tant to analyze the effects of imperfections and quan- 

■ turn errors on the algorithm accuracy. A usual model of 
quantum errors assumes that angles of unitary rotations 
fluctuates randomly in time for any qubit in some small 
interval e near the exact angle values determined by the 
ideal algorithm. In this case a realistic quantum com- 
putation remains close to the ideal one up to a number 
of performed gates N g ~ 1/e 2 . For example, the fidelity 
/ of computation, defined as a square of scalar product 
of quantum wavefunctions of ideal and perturbed algo- 
rithms, remains close to unity if a number of performed 
gates is smaller than N g . This result has been estab- 
lished analytically and numerically in extensive studies 
of various quantum algorithms dj ■ 

Another source of quantum errors comes from inter- 
nal imperfections generated by residual static couplings 

, between qubits and one-qubit energy level shifts which 
fluctuate from one qubit to another but remain static in 
time. These static imperfections may lead to appearance 
of many-body quantum chaos, which modifies strongly 

, the hardware properties of realistic quantum computer 
0]. The effects of static imperfections on the accuracy 
of quantum computation have been investigated on the 
examples of quantum algorithms for the models of com- 
plex quantum dynamics (see e.g. pH and Refs. therein) 
and a universal law for fidelity decay induced by static 
imperfections has been established [Hj for quantum al- 
gorithms simulating dynamics in the regime of quantum 
chaos. This law is based on the RMT treatment of imper- 
fections. At the same time it has been realized that the 
effects of static imperfections for dynamics in an inte- 
grable regime are not universal and more complicated. 
Therefore it is important to investigate the effects of 
static imperfections on an example of the well known 
Grover algorithm (GA) [S]. The results of these inves- 
tigations [Sj show that a quantum phase transition to 
quantum chaos takes place for the imperfection strength 
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Figure 1: Probability wc(t) of a searched state in GA 
as a function of the Grover iteration t at n to t = 12, e = 
0.002, to — 34.5. Left: curves show data for ideal GA, 
GA with gate to gate randomly fluctuating coefficients 
a.i, hj (see Eq. (2) in 0), GYQEC with l g = 10, GA with 
static imperfections (from top to bottom at t/t g = 1). 
Right: curves show data for GYQEC at l g — 1, 10, 20 and 
GA with static imperfections (also from top to bottom). 

e > e c ~ 1.7/ (rigy/ntot) where n g = O(n to t) is the to- 
tal number of quantum gates for one Grover iteration 
and ntot = n q + 1 is the total number of qubits. No- 
tations and detailed explanations are given in [S]. Here 
we give the description of the gyroscopic quantum error 
correction (GYQEC) method allowing significantly sup- 
press static imperfections in GA. For the first time this 
method is discussed in jj] for the quantum tent map. 

GYQEC is based on a random change of numeration 
of qubits after l g quantum gates. Namely, after l g gates 
about ntot I ^ swap operations are performed between ran- 
dom pairs of qubits so that the initial numeration of 
qubits is replaced by completely random one. However, 
in the quantum computer code this change is taken into 
account and the algorithm continues to run with new 
qubit numeration. In a sense the method uses a freedom 
of numeration of qubits in the program code and makes 
gyroscopic random rotations between all possibilities. 
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Figure 2: Husimi function in GA (compare with Fig. 2 
in 0), shown at moment t ~ to when wcit) has max- 
imum , for e = 0.002,0.004,0.008 (left to right respec- 
tively); n to t — 12. Top (bottom) row corresponds to the 
computation with (without) GYQEC at l g — 1. Den- 
sity is proportional to color changing from blue/zero to 
red/maximum. 




Figure 3: Left: probability wq, averaged over 10 disor- 
der realizations of static imperfections and taken at maxi- 
mum, for e = 0.002, 0.003, 0.004, 0.005, computation with 
GYQEC at different l g (top to bottom). Right: the gain 
factor R given by the ratio of wq (from left) to its max- 
imum value obtained in computations without GYQEC 
(same symbols). Here n to t = 12. 
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Figure 4: Probability Wi{t) in the four states (see Eq. (3) 
in |S]) as a function of iteration time t for n to t — 12, e = 
0.002, to — 34.5. From top to bottom, curves show data 
without GYQEC, with time fluctuating couplings, with 
GYQEC at l g = 1 (practically coincides with the previous 
curve) and at l g — 5, 10, 20. 



These rotations suppress the effects of static imperfec- 
tions (see Fig. 1). The GA accuracy is impoved with a 
decrease of l g even if at l g = 1 the GYQEC is not able to 
reach the case with randomly time fluctuating couplings 
between qubits JB|. A pictorial image of the accuracy 
improvement for the Husimi function in shown in Fig. 2. 
Indeed, GYQEC gives a significant increase of the proba- 
bility of searched state corresponding to lower horizontal 
line in a phase space square in Fig. 2. The variation of 
the searched probability wq with l g is shown in Fig. 3 for 
various values of e. GYQEC gives a maximal improve- 
ment of accuracy at minimal l g = 1 when the effect of 
randomization of static imperfections becomes maximal. 
For l g — 1 and n to t — 12 we reach the maximal accuracy 
gain factor R 6 which it is not very sensitive to e in 
a certain range. We expect that this R value will grow 
with the number of qubits n tot since in this case ran- 
dom rotations of computational basis will give stronger 
randomization of static imperfections. We note that the 
static imperfections preserve the total probability W4 in 
4-states (see 0) until e < e c while time fluctuations of 
couplings a,i,bij and GYQEC method give an exponen- 
tial time decay of W4 with a rate T oc e 2 (see Fig. 4). In 
spite of this decay we obtain the accuracy gain. 

In summary, we discussed here a new GYQEC method 
which performs random rotations in the computational 
basis of quantum computation code keeping track of these 
rotations in the quantum algorithm. This method uses 
a generic property of numeration freedom in the com- 
putational basis and allows to suppress significantly the 
effects of static imperfections. The GYQEC is rather 
general and can be applied to any quantum algorithm. 
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